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The  recent  synthesis  of  a  four-fold  silylated  sila-adamantane  molecule  (C24H72Sii4,  Td),  [J.  Fis¬ 
cher,  J.  Baumgartner,  and  C.  Marschner,  Science  310,  (2005)  825]  is  the  first  attempt  of  making 
the  silicon  analogue  of  adamantane.  It  has  adamantane-like  Siio  core,  capped  by  methyl  and  sily- 
methyl  ligands.  We  report  its  electronic  structure,  vibrational  spectrum,  and  the  infra-red  and 
Raman  spectra  calculated  within  the  density  functional  formalism  using  large  polarized  Gaussian 
basis  sets.  The  properties  of  sila-adamantane  are  compared  with  exact  silicon  analogue  (SiioHi4)  of 
adamantane.  Results  show  that  replacing  hydrogens  in  SiioHi6  by  methyl  and  silymethyl  ligands 
results  in  expansion  of  the  Siio  core  and  are  responsible  for  large  number  of  peaks  in  the  Raman 
spectrum.  The  Si-C  stretch  at  664  cm-1  and  methyl  deformations  frequencies  compare  well  with 
recent  measurements  of  vibrational  frequencies  of  methylated  silicon  surface. 
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I.  INTRODUCTION 


Silicon  is  of  considerable  interest  in  nanotechnology 
and  in  semiconductor  industry.  In  bulk  phase,  silicon 
has  diamond  lattice  like  carbon.  The  diamond  structure 
can  be  built  by  translating  the  smallest  unit  of  the  car¬ 
bon  lattice  called  adamantane.  The  adamantane  has  10 
carbon  atoms  that  are  capped  by  hydrogen  atoms  that 
passivate  the  dangling  bonds  of  carbon  atoms  in  the  sp3 
bonding.  The  adamantane  and  a  few  low-order  diamon- 
doids  have  been  synthesized  [ij.  Recently,  isolation  of 
several  new  lower-  and  medium-order  diamondoids  from 
petroleum  oil  was  reported  by  the  Molecular  Diamond 
Technologies  Group  at  Chevron-TexacoQ.  This  has  led 
to  a  renewal  of  interest  in  the  chemistry  and  physics  of 
the  diamondoids  Q. 

The  search  for  silicon  analogue  of  adamantane,  called 
sila-adamantane  SiioHig,  has  so  far  been  unsuccessful. 
There  has  been  however  some  progress.  A  recent  report 
by  Ficher,  Baumgartner  and  Marschner  (FBM)  of  syn¬ 
thesis  of  silylated  sila-adamantane  is  a  breakthrough  in 


the  synthesis  of  sila- adamantane  [4j.  Instead  of  hydrogen 
atoms  as  capping  elements  as  in  the  exact  silicon  ana¬ 
logue  of  adamantane,  silicon  atoms  in  sila-adamantane 
are  capped  by  methyl  and  trimethylsilyl  groups  (Cf.  Fig. 
[T]) .  It  was  characterized  using  the  nuclear  magnetic  reso¬ 
nance  (NMR)  spectroscopy.  The  ultra-violet  and  visible 
spectra  were  also  measured  which  showed  strong  absorp¬ 
tion  at  222  nm  (5.58  eV).  Followed  by  this  report,  Pic- 
cheri  reported  electronic  structure  of  C24Sii4H72  within 
Perdew-Burke-Ernzerhof  density  functional  model.  He 
reported  structural  parameters  calculated  using  the  3- 
21G*  basis  sets  and  density  fitting  methodology  [51- 
In  the  present  article,  we  report  the  electronic  struc¬ 
ture  of  sila-adamantane  examined  within  the  density 
functional  theory  using  large  polarized  Gaussian  basis 
sets.  Furthermore,  at  the  same  level  of  theory  we  cal¬ 
culate  vibrational  frequencies,  the  infra-red  and  Raman 
spectra.  We  also  calculate  these  properties  for  exact 
silicon  analogue  of  adamantane  (SiioH^),  and  compare 
them  with  those  of  sila-adamantane  C24Sii4H72  to  ex¬ 
amine  the  effect  of  methyl  and  silymethyl  ligands. 
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II.  FORMALISM 


Our  calculations  are  performed  within  the  Kohn-Sham 
formulation  jb]  of  the  density  functional  theory  Q  us¬ 
ing  NRLMOL  code  [9].  The  molecular  orbitals  are  ex¬ 
pressed  as  linear  combination  of  Gaussian  functions. 
The  exchange-correlation  effects  are  described  within  the 
generalized  gradient  approximation  using  Perdew-Burke- 
Ernzerhof[s]  parameterization.  For  a  reliable  prediction 
of  properties,  we  use  large  polarized  basis  sets  that  are 
specially  optimized  for  the  PBE  functional flp|.  Thus  the 
basis  for  Si  contains  6  s-,  5  p-,  and  3  d-type  Gaussians 
each  contracted  from  16  primitive  functions.  For  C,  12 
primitive  functions  were  used  to  contract  5  s-,  4p-,  and  3 
d-type  basis  functions.  The  basis  for  H  consist  of  4  s-,  3 
p-,  and  1  d-  type  basis  functions.  In  total  2614  basis  func¬ 
tions  were  used  in  calculation  of  C24Sii4H72-  Additional 
diffuse  functions  were  used  for  the  calculations  of  polar¬ 
izabilities.  The  exchange-correlation  contribution  to  ma¬ 
trix  elements  and  total  energy  are  obtained  numerically 
using  efficient  numerical  grids  employed  in  NRLMOL  [9(. 
The  vibrational  frequencies  were  obtained  by  diagonaliz¬ 
ing  the  dynamical  matrix  built  from  the  forces  computed 
at  perturbed  geometries.  Finally,  the  infra-red  and  Ra¬ 
man  spectra  were  computed  from  the  dipole  moment  and 
polarizability  derivatives,  respectively. 


FIG.  2:  (Color  online)  The  optimized  structure  of  C24Sii4H72 
without  hydrogens.  The  Si  atoms  that  are  part  of  the  sily- 
methyl  ligands  are  shown  in  dark  blue  for  clarity.  The  atomic 
labels  are  identical  to  those  in  experimental  structure  from 
Ref.  3. 


III.  RESULTS 

The  optimized  structure  of  SiioHi6  and  C24Sii4H72 
are  shown  in  Fig.  1.  The  C24Sii4H72  is  also  shown  with¬ 
out  hydrogens  in  Fig.  2.  To  distinguish  the  Si  atoms  that 
belong  to  silymethy  1-ligand,  these  four  atoms  are  shown 
in  dark  blue.  The  SiioHi6  core  in  C24Sii4H72  is  clearly 
seen.  The  atoms  have  labels  exactly  identical  to  those 
given  in  the  supplementary  information  of  the  experi¬ 
mental  study  in  Ref.  Q.  Both  structures  were  optimized 
within  the  T^  symmetry.  The  bond  distance  between  two 
silicon  atoms  (Si(l)-Si(2))  that  belong  to  the  tetrahedral 
adamantane  core  (Si(10))  is  2.39  A,  in  good  agreement 
with  experimental  value  of  2.36  A.  In  case  of  SiioHig, 
this  bond  length  is  calculated  to  be  2.36  A.  The  Si-SiCH3 
bond  length  (Si(2)-Si(5)  distance)  is  also  2.39  A  while  the 
experimental  value  is  2.35  A.  The  bond  angle  between 
the  Si-Si  bonds  from  the  Siio  core  (Si(2)-Si(l)-Si(2A))  in 
C24Sii4H72  is  110°  and  is  in  excellent  agreement  with  the 
experimental  value  of  110.3°.  Similarly,  very  good  agree- 
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ment  is  observed  for  the  dihedral  angle  between  three 
bonds  from  the  Siio  core  (  Si(2A)-Si(l)-Si(2)-Si(3)).  The 
predicted  angle  is  59.8°  while  experimental  value  is  59.9°. 
These  bond  angles  are  essentially  similar  in  SiioHi6  (an¬ 
gles  110.1°  and  59.  7°).  As  the  bond  angles  and  dihedral 
angles  in  SiioHi6  are  essentially  similar  to  those  in  the 
Siio  core  in  C24Sii4H72,  it  appears  that  replacing  hy¬ 
drogens  by  methyl  and  trimethylsilyl  ligands  leads  to  an 
expansion  of  the  Siio  core. 

IV.  ELECTRONIC  STRUCTURE  PROPERTIES 


The  fully  occupied  highest  occupied  molecular  or¬ 
bital  (HOMO)  of  C24Sii4H72  is  three  fold  degenerate 
and  has  T2  symmetry  while  its  LUMO  has  Ai  sym¬ 
metry.  The  HOMO-LUMO  (HL)  gap  is  4.16  eV.  For 
the  SiioHi6,  the  HL  gap  is  4.82  eV.  The  HL  gaps  re¬ 
ported  earlier  with  a  smaller  3-21G*  orbital  basis  and 
the  density  fitting  methodology  are  in  the  range  4.38- 
4.53  eV  depending  on  the  choice  of  exchange-correlation 
parameterization  [5].  We  also  calculated  the  HOMO- 
LUMO  gap  for  the  SiioHi6  using  the  time  dependent 
DFT  (TDDFT)  and  using  the  Perdew-Wang  GGA  as 


mg  t 
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implemented  in  ABINIT[ll].  The  calculated  gap  using 
TDDFT  of  4.6  eV  is  only  marginally  larger  than  the 
PBE-DFT  gap.  The  optical  gaps  of  SiioHi6  have  been 
studied  using  more  complex  quantum  Monte  Carlo  and 
GW  approximations  12  .  These  methods  give  value  of  9.2 
eV  for  the  optical  gaps.  Our  attempts  to  compute  the 
TDDFT  optical  gap  for  C24Sii4H72  were  not  successful 
due  to  large  memory  requirement. 


V.  VIBRATIONAL  SPECTRA 

A.  SiioHi6 

The  SiioHi6  cage  has  Td  point  group  symmetry.  It  has 
in  total  72  vibrational  modes  which  could  be  classified  as 
5Ai  +  A2  +  6  E  +  7Ti  +  11  T2.  Of  the  total  11  T2 
modes  that  are  infra-red  active,  only  four  have  apprecia¬ 
ble  intensity.  The  calculated  IR  spectra  is  shown  in  Fig. 
[3j  The  first  absorption  is  at  716  cm-1  (7.61  D/A2)  due 


FIG.  3:  (Color  online)  The  infra-red  spectra  of  SiioHi6  and 
C24Sil4H72. 


FIG.  4:  (Color  online)  The  Raman  spectra  of  SiioHi6. 


FIG.  5:  (Color  online)  The  Raman  spectra  of  C24Sii4H72. 
It  is  split  into  low  frequency  region  50-1000  cm-1  and  high 
frequency  2800-3100  cm-1  (CH  stretches)  region. 
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to  the  SiH2  wagging.  The  mode  at  885  cm-1  (3.85  D/A2) 
corresponds  to  the  bending  (scissor)  mode  of  SiH2.  The 
peaks  for  the  two  SiH  stretching  modes  are  at  higher  fre¬ 
quencies  of  2130  cm-1  and  2141  cm-1  are  coalesced  in 
the  figure.  The  quantities  in  bracket  are  intensities. 

The  Raman  spectra  of  SiioHi6  is  shown  in  Fig.  |4] 
There  are  5Ai,  6  E,  and  11  T2  Raman  active  modes. 
The  fully  symmetric  Ai  mode  at  2137  cm-1  is  the  sym¬ 
metric  SiH  stretch  and  is  the  most  intense  (intensity  130 
A4/amu).  The  absorption  due  to  fully  symmetric  stretch 
of  Si-Si  bonds  (breathing  of  silicon  cage)  occurs  at  lower 
(due  to  larger  mass)  frequency  of  299  cm-1  .  The  second 
most  intense  absorption  is  seen  at  lower  wave  number  of 
102  cm-1  with  intensity  of  68  A4/amu.  This  is  three  fold 
degenerate  Ti  mode  due  to  rocking  motion  of  SiH2.  The 
Raman  active  doubly  degenerate  mode  at  103  cm-1  is 
due  to  the  deformation  of  silicon  cage.  The  absorptions 
at  larger  frequencies  (  >  600  cm-1  )  involve  SiH2  defor¬ 
mations  and  SiH  stretches.  Less  intense  absorption  are 
620  cm-1  (SiH2  twist,  intensity  9  A4/amu),  688  cm-1 
(SiH2  and  SiH  rocking,  intensity  2  A4/amu),  and  878 
cm-1  (SiH2  scissor,  10  A4/amu  ). 

B.  C  24SL4H72 

We  computed  the  vibrational  frequencies,  within  the 
harmonic  approximation,  with  different  possible  orienta¬ 
tions  of  methyl  groups  that  are  consistent  with  the  Td 
symmetry.  In  all  cases  we  find  9  imaginary  frequen¬ 
cies  for  C24Sii4H72  due  to  the  rotation  of  CH3  group. 
These  rotations  seen  in  the  gas  phase  will  be  frozen  in 
the  crystalline  state.  The  total  324  vibrational  modes  in 
the  C24Sii4H72  can  be  classified  as  17A4  +  10  A2  +  27 
E  +  37T4  +  44  T2.  There  are  a  total  of  44  IR  active 
modes.  The  calculated  infra-red  absorption  spectrum  is 
displayed  in  Fig.  [3]  The  two  conspicuous  peaks  at  lower 
frequencies  of  795  cm-1  (8  D/A2)  and  834  cm-1  (  9 
D/A2)  are  due  to  the  SiCH3  rocking/wagging  motions. 
These  modes  are  similar  to  the  SiH2  wagging  and  rock¬ 
ing  seen  in  SiioHigThe  weak  absorption  at  1232  cm'1 
(  1.4  D/A2)  is  the  Si-CH3  umbrella  mode.  The  three 


peaks  at  higher  frequencies  in  the  frequency  range  2956- 
3032  cm-1  with  intensities  of  1.3-1. 4  D/A2  are  due  to 
the  symmetric  and  asymmetric  C-H  stretches.  A  number 
of  vibrational  modes  results  in  weak  absorption,  with  in¬ 
tensities  less  than  1  D/A2.  These  are  Si-C  asymmetric 
stretch  at  664  cm-1  (0.6  D/A2),  Si-Si-Si  deformation 
at  631  cm-1  (0.7  D/A2),  Si-C  symmetric  stretch  at  603 
cm-1  (0.6  D/A2),  CH3  deformation  at  1417  cm-1  (0.7 
D/A2),  and  C-H  symmetric  stretch  at  3060  (0.8  D/A2). 
The  compilation  of  the  frequencies  is  given  in  Table  HI 

TABLE  I:  The  IR  active  frequenters  for  the  C24Sii4H72.  The 
stretch  of  absorptions  is  given  in  bracket  with  following  ab¬ 
breviations:  w:  weak,  vs:  very  strong,  m:  medium.  The 
experimental  numbers  are  for  methylated  silicon(lll)  sur¬ 
face  obtained  by  high-resolution  electron  energy  loss  spec¬ 
troscopy  (HREELS)[3  an<i  transmission  infra-red  (TIR)Q 
spectroscopy. 


Freq.  (  cm  1  ' 

)  HREELS 

TIR  Mode  description 

664  (w) 

678 

Si-C  asymmetric  stretch 

603  (w) 

Si-C  symmetric  stretch 

631  (w) 

Si-Si-Si  deformation 

796(vs) 

SiCH3  bending 

834(vs) 

CH3  rocking 

1233(m) 

1257  CH3  Umbrella 

1417  (w) 

1424 

2856  CH3  deformation 

2956  (m) 

2909  C-H  symmetric  stretch 

2969  (m) 

2916 

2909  C-H  symmetric  stretch 

3032  (m) 

2989 

2928  C-H  asymmetric  stretch 

The  Raman  spectrum  of  C24Sii4H72  is  shown  in  Fig. 
[5]  The  spectrum  can  be  split  into  two  parts:  the  ab¬ 
sorptions  in  the  low-intermediate  frequency  range  (60- 
200  cm-1  )  and  intense  absorptions  in  the  frequency 
range  2900-3100  cm-1  .  The  doubly  degenerate  absorp¬ 
tion  at  59  cm-1  (77  A4/amu)  involves  rocking  like  mo¬ 
tion  of  Si-Si(CH3)3.  This  mode  is  depicted  in  Fig.  [6] 
The  absorption  at  84  cm-1  (63  A4/amu)  is  due  to  a 
doubly  degenerate  mode  that  results  in  deformation  of 
Si-Si-Si  bonds  in  the  Siio  core  similar  to  the  mode  at  103 
cm-1  seen  in  SiioHmThe  triplet  degenerate  mode  also 
due  to  the  distortion  of  Siio  core  at  85  cm-1  is  merged 
with  the  peak  at  84  cm-1  .  These  two  modes  are  shown 
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in  Fig.  □  The  weak  absorptions  at  108,  155  cm  1 2 *  are 
due  to  methyl  rotations.  The  most  intense  peak  in  the 
low  frequency  region  at  130  cm-1  (85  A4/amu)  is  due 
to  non-degenerate  Ai  mode  that  results  in  breathing  mo¬ 
tion  of  Si  atoms  to  which  Si (0113)3  ligands  are  attached. 
The  breathing  motions  of  rest  of  Si  atoms  in  Siio  core  (to 
which  two  CH3  groups  are  attached)  leads  to  an  absorp¬ 
tion  at  177  cm-1  (  A4/amu). 

The  intense  absorptions  at  higher  frequencies  are  due 
to  CH  stretches.  The  first  peak  at  2956  cm-1  (160 
A4/amu)  is  due  to  Ai  CH  stretch  where  CH3  belongs 
to  the  Si(CH3)3  group.  The  Ai  CH  stretch  in  the  two 
methyls  attached  to  silicon  atoms  leads  to  absorption  at 
slightly  higher  frequency  2969  cm-1  .  The  three- fold  de¬ 
generate  asymmetric  CH  stretches  absorb  at  3032-3033 
cm-1  (  102  A4/amu)  and  is  also  infra-red  active. 

The  most  obvious  difference  between  the  Raman  spec¬ 
tra  of  SiioHi6  and  C24Sii4H72  is  the  presence  of  the 
intense  peak  at  2130  cm-1  for  the  SiioHi6  involving 
SiH  bonds  and  the  high  frequency  modes  involving  CH 
bonds  in  the  methylated  C24Sii4H72  cluster.  It  may  be 
mentioned  here  that  there  are  several  very  low  inten¬ 
sity  modes  between  1000-1500  cm-1  for  the  C24Sii4H72. 
However,  the  intensities  of  these  peaks  are  comparatively 
very  small  and  hence  may  not  be  used  for  identification 
purposes. 

The  passivation  of  dangling  bonds  at  the  silicon  sur¬ 
faces  by  hydrogenation,  chlorination  and  alklation  (ethy¬ 
lation,  methylation  etc)  has  been  a  subject  of  sev¬ 
eral  studies  due  to  technological  importance  of  silicon. 
Two  groups  have  reported  experimental  studies  of  vi¬ 
brational  modes  of  methylated  silicon  surface  have  been 
reported  14] .  Although  the  sila-adamantane  core  studied 
here  is  a  very  small  unit  of  silicon  crystal,  it  is  inter¬ 


esting  to  compare  the  frequencies  due  to  Si-C  stretching 
and  methyl  deformations.  Yamada  and  coworkers  13  as¬ 
sign  the  vibration  of  678  cm-1  to  Si-C  which  compare 
well  with  664  cm-1  Si-C  asymmetric  stretch.  The  CH3 
stretch  and  deformation  freuenqcies  compare  well  (Cf.  |T|). 

VI.  CONCLUSION 


The  equilibrium  structure  of  recently  synthesized 
of  silylated  sila-adamantane  C24Sii4H72  was  obtained 
within  the  density  functional  theory,  using  large  polarized 
Gaussian  basis  sets  and  the  semi-local  PBE  approxima¬ 
tion  to  the  exchange-correlation  functional.  The  calcu¬ 
lated  geometrical  parameters  are  in  very  god  agreement 
with  experimental  values.  The  Siio  core  in  the  optimized 
C24Sii4H72  is  dilated  with  respected  to  the  silicon  ana¬ 
logue  of  adamantane  SiioHi6-  The  calculation  of  vibra¬ 
tional  frequencies  indicate  that  some  of  the  methyl  lig¬ 
ands  are  freely  rotating  in  the  gas  phase.  The  infra-red 
and  Raman  spectra  of  C24Sii4H72  is  reported  for  the  first 
time.  The  comparison  of  spectra  for  the  C24Sii4H72  with 
that  of  SiioHi6  do  not  show  any  similarity  since  in  both 
the  structures  only  the  Si  core  is  similar.  The  SiioHi6 
shows  intense  peaks  corresponding  to  stretching  of  the 
Si-H  bonds  which  are  missing  in  the  C24Sii4H72  due  to 
absence  of  SiH  bonds.  On  the  other  hand,  the  presence 
of  methyls  are  evident  in  the  low  frequency  peaks  for  the 
C24Sii4H72  and  also  in  very  high  frequency  peak  arising 
from  the  CH  bonds. 
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